Altered somatic hypermutation patterns in COVID-19 patients classifies disease severity.

Introduction: The success of the human body in fighting SARS-CoV2 infection relies on lymphocytes and their antigen receptors. Identifying and characterizing clinically relevant receptors is of utmost importance. Methods: We report here the application of a machine learning approach, utilizing B cell receptor repertoire sequencing data from severely and mildly infected individuals with SARS-CoV2 compared with uninfected controls. Results: In contrast to previous studies, our approach successfully stratifies non-infected from infected individuals, as well as disease level of severity. The features that drive this classification are based on somatic hypermutation patterns, and point to alterations in the somatic hypermutation process in COVID-19 patients. Discussion: These features may be used to build and adapt therapeutic strategies to COVID-19, in particular to quantitatively assess potential diagnostic and therapeutic antibodies. These results constitute a proof of concept for future epidemiological challenges.

B-Cell Receptor Features and Database Establishment in Recovered COVID-19 Patients by Combining 5'-RACE with PacBio Sequencing.

BACKGROUND: Antibodies induced by viral infection can not only prevent subsequent virus infection, but can also mediate pathological injury following infection. Therefore, understanding the B-cell receptor (BCR) repertoire of either specific neutralizing or pathological antibodies from patients convalescing from Coronavirus disease 2019 (COVID-19) infection is of benefit for the preparation of therapeutic or preventive antibodies, and may provide insight into the mechanisms of COVID-19 pathological injury. METHODS: In this study, we used a molecular approach of combining 5' Rapid Amplification of cDNA Ends (5'-RACE) with PacBio sequencing to analyze the BCR repertoire of all 5 IgH and 2 IgL genes in B-cells harvested from 35 convalescent patients after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. RESULTS: We observed numerous BCR clonotypes within most COVID-19 patients, but not in healthy controls, which validates the association of the disease with a prototypical immune response. In addition, many clonotypes were found to be frequently shared between different patients or different classes of antibodies. CONCLUSIONS: These convergent clonotypes provide a resource to identify potential therapeutic/prophylactic antibodies, or identify antibodies associated with pathological effects following infection with SARS-CoV-2.

The major role of junctional diversity in the horse antibody repertoire.

The antibody repertoire (Rep-seq) sequencing revolutionized the diversity of antigen B cell receptor studies, allowing deep and quantitative analysis to decipher the role of adaptive immunity in health and disease. Particularly, horse (Equus caballus) polyclonal antibodies have been produced and used since the century XIX to treat and prophylaxis diphtheria, tuberculosis, tetanus, pneumonia, and, more recently, COVID-19. However, our knowledge about the horse B cell receptors repertories is minimal. We present a deep horse antibody heavy chain repertoire (IGH) characterization of non-infected horses using NGS (Next generation sequencing). This study obtained a mean of 248,169 unique IgM clones and 66,141 unique IgG clones from four domestic adult horses. Rarefaction analysis showed sequence coverage was between 52 % and 82 % in IgM and IgG isotypes. We observed that besides horses antibody can use all functional IGHV genes, around 80 % of their antibodies use only three IGHV gene segments, and around 55 % use only one IGHJ gene segment. This limited VJ diversity seems to be compensated by the junctional diversity of these antibodies. We observed that the junctional diversity in horse antibodies is widespread, present in more than 90 % of horse antibodies. Besides this, the length of this region seems to be higher in horse antibodies than in other species. N1 and N2 nucleotides addition range from 0 to 111 nucleotides. In addition, around 45 % of the antibody clones have more than ten nucleotides in both the N1 and N2 junction regions. This diversity mechanism may be one of the most important in providing variability to the equine antibody repertoire. This study provides new insights regarding horse antibody composition, diversity generation, and particularities compared to other species, such as the frequency and length of N nucleotide addition. This study also points out the urgent need to better characterize TdT in horses and other species to better understand antibody repertoire characteristics.

AIRRscape: An interactive tool for exploring B-cell receptor repertoires and antibody responses.

The sequencing of antibody repertoires of B-cells at increasing coverage and depth has led to the identification of vast numbers of immunoglobulin heavy and light chains. However, the size and complexity of these Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) datasets makes it difficult to perform exploratory analyses. To aid in data exploration, we have developed AIRRscape, an R Shiny-based interactive web browser application that enables B-cell receptor (BCR) and antibody feature discovery through comparisons among multiple repertoires. Using AIRR-seq data as input, AIRRscape starts by aggregating and sorting repertoires into interactive and explorable bins of germline V-gene, germline J-gene, and CDR3 length, providing a high-level view of the entire repertoire. Interesting subsets of repertoires can be quickly identified and selected, and then network topologies of CDR3 motifs can be generated for further exploration. Here we demonstrate AIRRscape using patient BCR repertoires and sequences of published monoclonal antibodies to investigate patterns of humoral immunity to three viral pathogens: SARS-CoV-2, HIV-1, and DENV (dengue virus). AIRRscape reveals convergent antibody sequences among datasets for all three pathogens, although HIV-1 antibody datasets display limited convergence and idiosyncratic responses. We have made AIRRscape available as a web-based Shiny application, along with code on GitHub to encourage its open development and use by immuno-informaticians, virologists, immunologists, vaccine developers, and other scientists that are interested in exploring and comparing multiple immune receptor repertoires.

Novel skewed usage of B-cell receptors in COVID-19 patients with various clinical presentations.

B cell-mediated immune responses play important roles in controlling SARS-CoV infection. Here, we performed the single-cell B cell receptor sequencing (scBCR-seq) of the PBMC samples from eleven healthy controls, five asymptomatic subjects and 33 symptomatic COVID-19 patients with various clinical presentations, and subsequently analyzed the abundance and diversity of the BCR repertoires in different groups, respectively. We revealed the skewed usage of the IGHV, IGLV and IGKV genes and identified a number of heavy or light chain VDJ gene pairs and combinational preference in each group, such as IGKV3-7 and IGKV2-24 enriched in the asymptomatic subjects, whereas IGHV3-13, IGHV3-23-IGHJ4, IGHV1-18-IGLV3-19, IGHV1-18-IGLV3-21, and IGHV1-18-IGLV3-25 enriched in the recovery patients with severe diseases. We also observed the differential expression of IGHV3-23 in various B cell clusters by analysis of the scRNA-seq data. Additional dock analysis indicated that IGHV3-13 could bind to the spike protein of SARS-CoV-2. These findings may advance our understanding of the humoral immune responses in COVID-19 patients and help develop novel vaccine candidates as well as therapeutical antibodies against SASR-CoV-2 infections.

Establishment of the large-scale longitudinal multi-omics dataset in COVID-19 patients: data profile and biospecimen.

Understanding and monitoring virus-mediated infections has gained importance since the global outbreak of the coronavirus disease 2019 (COVID-19) pandemic. Studies of high-throughput omics-based immune profiling of COVID-19 patients can help manage the current pandemic and future virus-mediated pandemics. Although COVID-19 is being studied since past 2 years, detailed mechanisms of the initial induction of dynamic immune responses or the molecular mechanisms that characterize disease progression remains unclear. This study involved comprehensively collected biospecimens and longitudinal multi-omics data of 300 COVID-19 patients and 120 healthy controls, including whole genome sequencing (WGS), single-cell RNA sequencing combined with T cell receptor (TCR) and B cell receptor (BCR) sequencing (scRNA(+scTCR/BCR)-seq), bulk BCR and TCR sequencing (bulk TCR/BCR-seq), and cytokine profiling. Clinical data were also collected from hospitalized COVID-19 patients, and HLA typing, laboratory characteristics, and COVID-19 viral genome sequencing were performed during the initial diagnosis. The entire set of biospecimens and multi-omics data generated in this project can be accessed by researchers from the National Biobank of Korea with prior approval. This distribution of largescale multi-omics data of COVID-19 patients can facilitate the understanding of biological crosstalk involved in COVID-19 infection and contribute to the development of potential methodologies for its diagnosis and treatment. [BMB Reports 2022; 55(9): 465-471].

Comparative global B cell receptor repertoire difference induced by SARS-CoV-2 infection or vaccination via single-cell V(D)J sequencing.

Dynamic changes of the paired heavy and light chain B cell receptor (BCR) repertoire provide an essential insight into understanding the humoral immune response post-SARS-CoV-2 infection and vaccination. However, differences between the endogenous paired BCR repertoire kinetics in SARS-CoV-2 infection and previously recovered/naïve subjects treated with the inactivated vaccine remain largely unknown. We performed single-cell V(D)J sequencing of B cells from six healthy donors with three shots of inactivated SARS-CoV-2 vaccine (BBIBP-CorV), five people who received the BBIBP-CorV vaccine after having recovered from COVID-19, five unvaccinated COVID-19 recovered patients and then integrated with public data of B cells from four SARS-CoV-2-infected subjects. We discovered that BCR variable (V) genes were more prominently used in the SARS-CoV-2 exposed groups (both in the group with active infection and in the group that had recovered) than in the vaccinated groups. The VH gene that expanded the most after SARS-CoV-2 infection was IGHV3-33, while IGHV3-23 in the vaccinated groups. SARS-CoV-2-infected group enhanced more BCR clonal expansion and somatic hypermutation than the vaccinated healthy group. A small proportion of public clonotypes were shared between the SARS-CoV-2 infected, vaccinated healthy, and recovered groups. Moreover, several public antibodies had been identified against SARS-CoV-2 spike protein. We comprehensively characterize the paired heavy and light chain BCR repertoire from SARS-CoV-2 infection to vaccination, providing further guidance for the development of the next-generation precision vaccine.

Optimization of peripheral blood volume for in silico reconstitution of the human B-cell receptor repertoire.

B cells recognize antigens via membrane-expressed B-cell receptors (BCR) and antibodies. Similar human BCR sequences are frequently found at a significantly higher frequency than that theoretically calculated. Patients infected with SARS-CoV2 and HIV or with autoimmune diseases share very similar BCRs. Therefore, in silico reconstitution of BCR repertoires and identification of stereotypical BCR sequences related to human pathology have diagnostic potential. Furthermore, monitoring changes of clinically significant BCR sequences and isotype conversion has prognostic potential. For BCR repertoire analysis, peripheral blood (PB) is the most convenient source. However, the optimal human PB volume for in silico reconstitution of the BCR repertoire has not been studied in detail. Here, we sampled 5, 10, and 20 mL PB from the left arm and 40 mL PB from the right arm of two volunteers, reconstituted in silico PB BCR repertoires, and compared their composition. In both volunteers, PB sampling over 20 mL resulted in slight increases in functional unique sequences (FUSs) or almost no increase in repertoire diversity. All FUSs with a frequency above 0.08% or 0.03% in the 40 mL PB BCR repertoire were detected even in the 5 mL PB BCR repertoire from each volunteer. FUSs with a higher frequency were more likely to be found in BCR repertoires from reduced PB volume, and those coexisting in two repertoires showed a statistically significant correlation in frequency irrespective of sampled anatomical site. The correlation was more significant in higher-frequency FUSs. These observations support the potential of BCR repertoire analysis for diagnosis.

Landscape of B Cell Receptor Repertoires in COVID-19 Patients Revealed Through CDR3 Sequencing of Immunoglobulin Heavy and Light Chains.

The outbreak and persistence of coronavirus disease 2019 (COVID-19) threaten human health. B cells play a vital role in fighting the infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Despite many studies on the immune responses in COVID-19 patients, it is still unclear how B cell receptor (BCR) constituents, including immunoglobulin heavy (IGHs) and light chains (IGLs), respond to SARS-CoV-2 in patients with varying symptoms. In this study, we conducted complementarity-determining region 3 (CDR3) sequencing of BCR IGHs and IGLs from the peripheral blood of COVID-19 patients and healthy donors. The results showed significantly reduced clonal diversity, more expanded clones, and longer CDR3 lengths of IGH and IGL in COVID-19 patients than those in healthy individuals. The IGLs had a much higher percentage of VJ skew usage (47.83% IGLV and 42.86% IGLJ were significantly regulated) than the IGHs (12.09% IGHV and 0% IGHJ) between the healthy individuals and patients, which indicated the importance of BCR light chains. Furthermore, we found a largely expanded IGLV3-25 gene cluster mostly pairing with IGLJ1 and ILGJ2 in COVID-19 patients and a newly identified upregulated IGLJ1 gene and IGLJ2+IGLV13-21 recombination, both of which are potential sources of SARS-CoV-2-targeting antibodies. Our findings on specific immune B-cell signatures associated with COVID-19 have clinical implications for vaccine and biomarker development for disease diagnosis.

Analysis of B Cell Receptor Repertoires Reveals Key Signatures of the Systemic B Cell Response after SARS-CoV-2 Infection.

A comprehensive study of the B cell response against SARS-CoV-2 could be significant for understanding the immune response and developing therapeutical antibodies and vaccines. To define the dynamics and characteristics of the antibody repertoire following SARS-CoV-2 infection, we analyzed the mRNA transcripts of immunoglobulin heavy chain (IgH) repertoires of 24 peripheral blood samples collected between 3 and 111 days after symptom onset from 10 COVID-19 patients. Massive clonal expansion of naive B cells with limited somatic hypermutation (SHM) was observed in the second week after symptom onset. The proportion of low-SHM IgG clones strongly correlated with spike-specific IgG antibody titers, highlighting the significant activation of naive B cells in response to a novel virus infection. The antibody isotype switching landscape showed a transient IgA surge in the first week after symptom onset, followed by a sustained IgG elevation that lasted for at least 3 months. SARS-CoV-2 infection elicited poly-germ line reactive antibody responses. Interestingly, 17 different IGHV germ line genes recombined with IGHJ6 showed significant clonal expansion. By comparing the IgH repertoires that we sequenced with the 774 reported SARS-CoV-2-reactive monoclonal antibodies (MAbs), 13 shared spike-specific IgH clusters were found. These shared spike-specific IgH clusters are derived from the same lineage of several recently published neutralizing MAbs, including CC12.1, CC12.3, C102, REGN10977, and 4A8. Furthermore, identical spike-specific IgH sequences were found in different COVID-19 patients, suggesting a highly convergent antibody response to SARS-CoV-2. Our analysis based on sequencing antibody repertoires from different individuals revealed key signatures of the systemic B cell response induced by SARS-CoV-2 infection. IMPORTANCE Although the canonical delineation of serum antibody responses following SARS-CoV-2 infection has been well established, the dynamics of antibody repertoire at the mRNA transcriptional level has not been well understood, especially the correlation between serum antibody titers and the antibody mRNA transcripts. In this study, we analyzed the IgH transcripts and characterized the B cell clonal expansion and differentiation, isotype switching, and somatic hypermutation in COVID-19 patients. This study provided insights at the repertoire level for the B cell response after SARS-CoV-2 infection.

Efficient discovery of SARS-CoV-2-neutralizing antibodies via B cell receptor sequencing and ligand blocking.

Although several monoclonal antibodies (mAbs) targeting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been approved for coronavirus disease 2019 (COVID-19) therapy, development was generally inefficient, with lead generation often requiring the production and testing of numerous antibody candidates. Here, we report that the integration of target-ligand blocking with a previously described B cell receptor-sequencing approach (linking B cell receptor to antigen specificity through sequencing (LIBRA-seq)) enables the rapid and efficient identification of multiple neutralizing mAbs that prevent the binding of SARS-CoV-2 spike (S) protein to angiotensin-converting enzyme 2 (ACE2). The combination of target-ligand blocking and high-throughput antibody sequencing promises to increase the throughput of programs aimed at discovering new neutralizing antibodies.

B cell receptor repertoire kinetics after SARS-CoV-2 infection and vaccination.

B cells are important in immunity to both severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection and vaccination, but B cell receptor (BCR) repertoire development in these contexts has not been compared. We analyze serial samples from 171 SARS-CoV-2-infected individuals and 63 vaccine recipients and find the global BCR repertoire differs between them. Following infection, immunoglobulin (Ig)G1/3 and IgA1 BCRs increase, somatic hypermutation (SHM) decreases, and, in severe disease, IgM and IgA clones are expanded. In contrast, after vaccination, the proportion of IgD/M BCRs increase, SHM is unchanged, and expansion of IgG clones is prominent. VH1-24, which targets the N-terminal domain (NTD) and contributes to neutralization, is expanded post infection except in the most severe disease. Infection generates a broad distribution of SARS-CoV-2-specific clones predicted to target the spike protein, while a more focused response after vaccination mainly targets the spike's receptor-binding domain. Thus, the nature of SARS-CoV-2 exposure differentially affects BCR repertoire development, potentially informing vaccine strategies.

Single-cell multi-omics reveals dyssynchrony of the innate and adaptive immune system in progressive COVID-19.

Dysregulated immune responses against the SARS-CoV-2 virus are instrumental in severe COVID-19. However, the immune signatures associated with immunopathology are poorly understood. Here we use multi-omics single-cell analysis to probe the dynamic immune responses in hospitalized patients with stable or progressive course of COVID-19, explore V(D)J repertoires, and assess the cellular effects of tocilizumab. Coordinated profiling of gene expression and cell lineage protein markers shows that S100Ahi/HLA-DRlo classical monocytes and activated LAG-3hi T cells are hallmarks of progressive disease and highlights the abnormal MHC-II/LAG-3 interaction on myeloid and T cells, respectively. We also find skewed T cell receptor repertories in expanded effector CD8+ clones, unmutated IGHG+ B cell clones, and mutated B cell clones with stable somatic hypermutation frequency over time. In conclusion, our in-depth immune profiling reveals dyssynchrony of the innate and adaptive immune interaction in progressive COVID-19.

The autoimmune signature of hyperinflammatory multisystem inflammatory syndrome in children.

Multisystem inflammatory syndrome in children (MIS-C) manifests as a severe and uncontrolled inflammatory response with multiorgan involvement, occurring weeks after SARS-CoV-2 infection. Here, we utilized proteomics, RNA sequencing, autoantibody arrays, and B cell receptor (BCR) repertoire analysis to characterize MIS-C immunopathogenesis and identify factors contributing to severe manifestations and intensive care unit admission. Inflammation markers, humoral immune responses, neutrophil activation, and complement and coagulation pathways were highly enriched in MIS-C patient serum, with a more hyperinflammatory profile in severe than in mild MIS-C cases. We identified a strong autoimmune signature in MIS-C, with autoantibodies targeted to both ubiquitously expressed and tissue-specific antigens, suggesting autoantigen release and excessive antigenic drive may result from systemic tissue damage. We further identified a cluster of patients with enhanced neutrophil responses as well as high anti-Spike IgG and autoantibody titers. BCR sequencing of these patients identified a strong imprint of antigenic drive with substantial BCR sequence connectivity and usage of autoimmunity-associated immunoglobulin heavy chain variable region (IGHV) genes. This cluster was linked to a TRBV11-2 expanded T cell receptor (TCR) repertoire, consistent with previous studies indicating a superantigen-driven pathogenic process. Overall, we identify a combination of pathogenic pathways that culminate in MIS-C and may inform treatment.

SARS-CoV-2-Reactive Mucosal B Cells in the Upper Respiratory Tract of Uninfected Individuals.

SARS-CoV-2 is a respiratory pathogen that can cause severe disease in at-risk populations but results in asymptomatic infections or a mild course of disease in the majority of cases. We report the identification of SARS-CoV-2-reactive B cells in human tonsillar tissue obtained from children who were negative for coronavirus disease 2019 prior to the pandemic and the generation of mAbs recognizing the SARS-CoV-2 Spike protein from these B cells. These Abs showed reduced binding to Spike proteins of SARS-CoV-2 variants and did not recognize Spike proteins of endemic coronaviruses, but subsets reacted with commensal microbiota and exhibited SARS-CoV-2-neutralizing potential. Our study demonstrates pre-existing SARS-CoV-2-reactive Abs in various B cell populations in the upper respiratory tract lymphoid tissue that may lead to the rapid engagement of the pathogen and contribute to prevent manifestations of symptomatic or severe disease.

Single-cell immune profiling reveals distinct immune response in asymptomatic COVID-19 patients.

While some individuals infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) present mild-to-severe disease, many SARS-CoV-2-infected individuals are asymptomatic. We sought to identify the distinction of immune response between asymptomatic and moderate patients. We performed single-cell transcriptome and T-cell/B-cell receptor (TCR/BCR) sequencing in 37 longitudinal collected peripheral blood mononuclear cell samples from asymptomatic, moderate, and severe patients with healthy controls. Asymptomatic patients displayed increased CD56briCD16- natural killer (NK) cells and upregulation of interferon-gamma in effector CD4+ and CD8+ T cells and NK cells. They showed more robust TCR clonal expansion, especially in effector CD4+ T cells, but lack strong BCR clonal expansion compared to moderate patients. Moreover, asymptomatic patients have lower interferon-stimulated genes (ISGs) expression in general but large interpatient variability, whereas moderate patients showed various magnitude and temporal dynamics of the ISGs expression across multiple cell populations but lower than a patient with severe disease. Our data provide evidence of different immune signatures to SARS-CoV-2 in asymptomatic infections.

A Single-Cell Atlas of Lymphocyte Adaptive Immune Repertoires and Transcriptomes Reveals Age-Related Differences in Convalescent COVID-19 Patients.

COVID-19 disease outcome is highly dependent on adaptive immunity from T and B lymphocytes, which play a critical role in the control, clearance and long-term protection against SARS-CoV-2. To date, there is limited knowledge on the composition of the T and B cell immune receptor repertoires [T cell receptors (TCRs) and B cell receptors (BCRs)] and transcriptomes in convalescent COVID-19 patients of different age groups. Here, we utilize single-cell sequencing (scSeq) of lymphocyte immune repertoires and transcriptomes to quantitatively profile the adaptive immune response in COVID-19 patients of varying age. We discovered highly expanded T and B cells in multiple patients, with the most expanded clonotypes coming from the effector CD8+ T cell population. Highly expanded CD8+ and CD4+ T cell clones show elevated markers of cytotoxicity (CD8: PRF1, GZMH, GNLY; CD4: GZMA), whereas clonally expanded B cells show markers of transition into the plasma cell state and activation across patients. By comparing young and old convalescent COVID-19 patients (mean ages = 31 and 66.8 years, respectively), we found that clonally expanded B cells in young patients were predominantly of the IgA isotype and their BCRs had incurred higher levels of somatic hypermutation than elderly patients. In conclusion, our scSeq analysis defines the adaptive immune repertoire and transcriptome in convalescent COVID-19 patients and shows important age-related differences implicated in immunity against SARS-CoV-2.

Global characterization of B cell receptor repertoire in COVID-19 patients by single-cell V(D)J sequencing.

The world is facing a pandemic of Corona Virus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Adaptive immune responses are essential for SARS-CoV-2 virus clearance. Although a large body of studies have been conducted to investigate the immune mechanism in COVID-19 patients, we still lack a comprehensive understanding of the BCR repertoire in patients. In this study, we used the single-cell V(D)J sequencing to characterize the BCR repertoire across convalescent COVID-19 patients. We observed that the BCR diversity was significantly reduced in disease compared with healthy controls. And BCRs tend to skew toward different V gene segments in COVID-19 and healthy controls. The CDR3 sequences of heavy chain in clonal BCRs in patients were more convergent than that in healthy controls. In addition, we discovered increased IgG and IgA isotypes in the disease, including IgG1, IgG3 and IgA1. In all clonal BCRs, IgG isotypes had the most frequent class switch recombination events and the highest somatic hypermutation rate, especially IgG3. Moreover, we found that an IgG3 cluster from different clonal groups had the same IGHV, IGHJ and CDR3 sequences (IGHV4-4-CARLANTNQFYDSSSYLNAMDVW-IGHJ6). Overall, our study provides a comprehensive characterization of the BCR repertoire in COVID-19 patients, which contributes to the understanding of the mechanism for the immune response to SARS-CoV-2 infection.

Dynamics of B cell repertoires and emergence of cross-reactive responses in patients with different severities of COVID-19.

Individuals with the 2019 coronavirus disease (COVID-19) show varying severity of the disease, ranging from asymptomatic to requiring intensive care. Although monoclonal antibodies specific to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been identified, we still lack an understanding of the overall landscape of B cell receptor (BCR) repertoires in individuals with COVID-19. We use high-throughput sequencing of bulk and plasma B cells collected at multiple time points during infection to characterize signatures of the B cell response to SARS-CoV-2 in 19 individuals. Using principled statistical approaches, we associate differential features of BCRs with different disease severity. We identify 38 significantly expanded clonal lineages shared among individuals as candidates for responses specific to SARS-CoV-2. Using single-cell sequencing, we verify the reactivity of BCRs shared among individuals to SARS-CoV-2 epitopes. Moreover, we identify the natural emergence of a BCR with cross-reactivity to SARS-CoV-1 and SARS-CoV-2 in some individuals. Our results provide insights important for development of rational therapies and vaccines against COVID-19.

Germline IGHV3-53-encoded RBD-targeting neutralizing antibodies are commonly present in the antibody repertoires of COVID-19 patients.

Monoclonal antibodies (mAbs) encoded by IGHV3-53 (VH3-53) targeting the spike receptor-binding domain (RBD) have been isolated from different COVID-19 patients. However, the existence and prevalence of shared VH3-53-encoded antibodies in the antibody repertoires is not clear. Using antibody repertoire sequencing, we found that the usage of VH3-53 increased after SARS-CoV-2 infection. A highly shared VH3-53-J6 clonotype was identified in 9 out of 13 COVID-19 patients. This clonotype was derived from convergent gene rearrangements with few somatic hypermutations and was evolutionary conserved. We synthesized 34 repertoire-deduced novel VH3-53-J6 heavy chains and paired with a common IGKV1-9 light chain to produce recombinant mAbs. Most of these recombinant mAbs (23/34) possess RBD binding and virus-neutralizing activities, and recognize ACE2 binding site via the same molecular interface. Our computational analysis, validated by laboratory experiments, revealed that VH3-53 antibodies targeting RBD are commonly present in COVID-19 patients' antibody repertoires, indicating many people have germline-like precursor sequences to rapidly generate SARS-CoV-2 neutralizing antibodies. Moreover, antigen-specific mAbs can be digitally obtained through antibody repertoire sequencing and computational analysis.